%% FigS21A
load('FigS21\G0012K - Full_A');
load('FigS21\G0012K - Half_A');
load('FigS21\G0012K - Halfn_A');
Vmf = data12fl(:,2);
Vsf = data12fl(:,1);
Vmhn = data12h(:,2);
Vshn = data12h(:,1);
Vmhp = data12hn(:,2);
Vshp = data12hn(:,1);

figure(10)
plot(Vsf(460:800)./(10*99.5e-5),Vmf(460:800).*10);
xlim([-100 100]);ylim([-25 25])
figure(11)
plot(Vsf(800:1140)./(10*99.5e-5),Vmf(800:1140).*10);
xlim([-100 100]);ylim([-25 25])
figure(12)
plot(Vshn(330:850)./(10*99.5e-5),Vmhn(330:850).*10);
xlim([-100 100]);ylim([-25 25])
figure(13)
plot(Vshp(2600:3250)./(10*99.5e-5),Vmhp(2600:3250).*10)
xlim([-100 100]);ylim([-25 25])

%% FigS21B
load(strcat('FigS21\C020 - squarewave 4_3K 400Hz'))
Vmf = fswp4K(:,2);
Vsf = fswp4K(:,1);
figure(120)
plot((1:14000)./270000,Vsf(1:14000)./(10*99.5e-5));

hold on; 
plot((1:14000)./270000,Vmf(1:14000)./(20*99.5e-5));
hold off
xlim([0 0.05])